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ABSTRACT 

The dynamical effects of magnetic fields in models of radiative, Herbig-Haro 
(HH) jets have been studied in a number of papers. For example, magnetized, 
radiative jets from variable sources have been studied with axisymmetric and 
3D numerical simulations. In this paper, we present an analytic model describ- 
ing the effect of a toroidal magnetic field on the internal working surfaces that 
result from a variability in the ejection velocity. We find that for parameters 
appropriate for HH jets the forces associated with the magnetic field dominate 
over the gas pressure force within the working surfaces. Depending on the ram 
pressure radial cross section of the jet, the magnetic field can produce a strong 
axial pinch, or, alternatively, a broadening of the internal working surfaces. We 
check the validity of the analytic model with axisymmetric numerical simulations 
of variable, magnetized jets. 

Subject headings: ISM: kinematics and dynamics - ISM: jets and outflows - ISM: 
Herbig-Haro objects - stars: magnetic fields - stars: pre-main sequence - winds, 
outflows 

1. Introduction 



It is now relatively certain that some Herbig-Haro (HH) jets have knot structures which 
are the result of a time-variability in the ejection. For example, the observations of some 
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jets with organized structures of kn ots of different sizes (e. g ., HH 30, 34 and 111, see 
Esquivel et al.l 120071 . iRaga et al.ll2002l and iMasciadri et al.ll2002l ) can be reproduced surpris- 
ingly well with variable ejection jet models. In the present paper, we study the effect of the 
presence of a magnetic field on the evolution of a variable jet. 

It is still an open question to what extent magnetic fields are important in determining 

the dynamics of HH jets. The associ ated problem of radiativ e, M HD jets has been explored in 

some detail in the existing literature. ICerqueira et al.l (119971 ). and lCerqueira fc de Gouveia Dal Pino 
(119991 ) computed 3D simulations of radiative, MHD l ets with different magnetic field config- 
urations (at the injection point). iFrank et al.l (119981 ) carried out axisymmetric simulations 
of similar flows. 

The problem of an MHD, radiative jet ejected with a time-v a riable velocity was ex - 
plored with axisymme t ric simulations by Gardiner fc Frankl (200oj) ; Gardiner et al. I dgQOo]); 
Stone fc Harde3 J2000h:lo , Sullivan fc Ravi J200oh : lFrank et al.l (l200oh : be Colle fc RagJj2006h 



and 



Hartigan et al.l (2007). Variable, MH D jets were also explored with 3D simulations by 



Cerqueira fc de Gouveia Dal Pino ( 2001al lbl). The general conclusions that can be obtained 
from these simulations is that the internal working surfaces produced by the ejection vari- 
ability are not affected strongly by a poloidal magnetic field. On the other hand, if the 
magnetic field is toroidal (or, alternatively, has a strong toroidal component), the material 
within the working surfaces of the jet flow has a stronger concentration towards the jet axis. 



Gardiner &: Frankl (120001 ) showed that in a variable ejection velocity jet the "continuous 



jet beam" sections in between the working surfaces have a low toroidal magnetic field, which 
grows in strength quite dramatically when the material goes through one of the working 
surface shocks into one of the knots. In the present paper, we present a simple, analytic 
model from which we obtain the conditions under which the toroidal magnetic field produces 
an axial compression of the internal working surfaces. This analytic model is presented in 
§2. In §3, we present axisymmetric numerical simulations in which we compare the working 
surfaces with and without a toroidal magnetic field, showing the effect described by the 
analytic model. Finally, in §4 we present our conclusions. 



2. The radial motion of the material within an internal working surface 



2.1. General considerations 



A time-variability in the ejection velocity leads to the formation of two-shock "internal 
working surfaces" which travel down the jet flow. In a frame of reference that moves with the 
working surface, the flow takes the configuration shown in Figure 1, with material entering 
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the shocked layer from both the upstream and downstream directions. 

Let us consider an internal working surface within a cylindrically symmetric jet with a 
toroidal magnetic field configuration. The material in the jet beam cross section or within 
the working surface is subject to two radial forces: the magnetic pinch force 

° d (rB), (I) 



4irr dr 



where B is the toroidal magnetic field and r the cylindrical radius, and the force due to the 
pressure gradient 

dP 

where P is the gas pressure. The cross section of the jet is in lateral equilibrium when 
F = F m + F p = 0, it will be subject to a lateral expansion when F > and to a compression 
when F < 0. 

Let us now assume that the jet beam has a generic cross section of the form 

P{r) = Pop(r) , (3) 

B(r) = B B(r ) , (4) 

v(r) = Vov(r) , (5) 

where p(r) is the density, B{r) the magnitude of the (toroidal) magnetic field, and v(r) = 
Vj—v w is the relative velocity with which the jet material (moving at a velocity Vj) enters the 
working surface (which moves at a velocity v w ), see Figure 1. The constants p , B Q and v 
correspond to characteristic values of the respective quantities, and p(r), B_{r) and v(r) are 
dimensionless functions of the radius giving the radial dependence of the flow variables from 
r = (the symmetry axis) out to r = rj (the outer radius of the jet beam). In principle, 
these three dimensionless functions are of order one unless very strong changes in the flow 
variables occur across the jet cross section. 

Let us now consider that the material goes through the "Mach disk" shock of an internal 
working surface. If we assume that the shock is strong (i. e., that it is highly supersonic 



and superalfvenic), from the Rankine-Hugoniot equations for MHD (e. g. iDraine McKee 



19931 ) the post-shock radial cross section is given by : 



P: r (r) = ^-p Q v 2 p(r)v 2 (r); 
7 + 1 - 



prad(r) _ (8^ 2 p 3 /2 v clP 3/2 (r)v(r) 
w K) B B{r) ' 1 ' 
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7 + 1 

Pl r ( r ) = tPopJV) I 

P ™ 1 J P S(r) ' 1 ' 

BZ r (r) = ^-B B(r); 
7-1 

c d w = (8^ i/2 p; /2 v /2 (o^), (8) 

where P™ r (r), pZ(r) and B™ r (r) are the post-shock gas pressure, density and magnetic field 
cross sections for the case of a non-radiative shock, and P™ (r), p^ ad (V) and B™ d (r) are the 
cross sections for the case of a radiative shock in which the post-shock gas instantaneously 
cools to an isothermal sound speed c w . As we have said above, equations ([ME]) have been 
derived for the case of a strong shock. In order to obtain these relations it is also necessary to 
assume that the pre-shock Alfvenic Mach number has values smaller than ~ M 2 = (v/c w ) 2 . 

The factors including the specific heat ratio 7 (see equations [MS]) take the numerical 
values 2/(7+1) = 3/4 and (7 + l)/(7 — 1) = 4 for the case of a monoatomic, non-relativistic 
gas (i. e., for 7 = 5/3). From now on, we will use these numerical values in order to simplify 
the equations. 

Combining equations (JMH1) with (P12]) we obtain the magnetic and gas pressure forces 
acting radially on the post-Mach disk material. The resulting magnetic force is 

K r = —f%(r) ; F™ d =^f™ d (r), (9) 
irrj 77 

where the dimensionless force f m (r) is given by 

/:« = -5(r)^[r5(r)] ; 

f£*(r) = -^(r) 2 (r)^ [rp 1/2 (r)v(r)] . (10) 



The resulting gas pressure force is 



3P0^| f nr M . F rad _ i^) 1 ^ 'ftj^ Upcg, , rad 

r, Ip [r) ' ^ r,P 
where the dimensionless force f p (r) is given by 



= -j^fZ r (r) ; P; ad = v 7 ^ "" /TV), (11) 



W = fe(0- 2 (0] ; f; ad (r) = -r± 



p 3 / 2 (r)v(r) 



B(r) 



(12) 



- 5 - 



2.2. Scaling properties of the magnetic and gas pressure forces 

Let us now consider the ratio Mj P between the moduli of the magnetic and gas pressure 
forces. From equations (Q and ffTTj) we obtain 



(M/P) nr 



64 
3Mj 



/7( 



(M/P) 



rad 



V2M- 
M A 



2 



fT(r) 



fl ad (r) 



(13) 



where Ma = vo/va is the Alfvenic Mach number (obtained with the characteristic velocity v$ 
and the Alfven velocity va = B / \/ 'Aix p ) , M w = v /c w is the sonic Mach number (calculated 
with the characteristic velocity v and the post-shock sound speed c w of the radiative shock) 
and the / m (r) and f p (r) functions are given by equations (fTOl and f|T2|) . respectively 

One can argue that if the dimensionless cross section of the jet (described by equations 
[3KHI) is smooth, then the f^ir), fm d ( r ) an d fp( r ) functions (see equations fTUI and fT2l) will 
have values of order 1. 

In our derivation of the pressure force within the internal working surface, we have only 
considered the gradient of the post-shock gas pressure. Of course, the fact that the working 
surface material is free to leave through the sides of the jet beam will lead to an extra gas 
pressure gradient (directed outwards), particularly in the case of a non-radiative flow. The 
dimensionless pressure cross section due to this effect is still likely to lead to a dimensionless 
force f p {r) ~ 1. 

Setting fmi r )i fm d ( r )i fp r { r )-> fp ad ( r ) ~ I; from equation (1T3|) we then obtain 



< M ^~3^r (M/PU ~fj- (14) 

From these two estimates of the ratio between the magnetic and gas pressure forces, we 
conclude that 



• for the non-radiative case : if the Alvenic Mach number of the flow entering the Mach 
disk is large (e. g., Ma > 10) we have (M/P) nr <C 1, and therefore the lateral expansion 
or contraction of the gas within the working surface will be governed by the gas pressure 
force, 

• for the radiative case : if we consider jets with given values for va and c w , it is clear 
that as the velocity v increases, the (M/P) rad ratio increases (proportional to v ). 
In particular, if we have flows with v A ~ c w , the magnetic to gas pressure force ratio 
has values (M/P) ra d ~ M w . Thus, for a Mach disk in the strong shock regime, the 
post-shock magnetic pressure force will under most conditions dominate over the gas 
pressure force. 
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Therefore, for the non-radiative and the radiative cases, whether the jet material within the 
working surface expands or contracts in the radial direction will be determined by the signs 
of fp T (r) and f™ d (r), respectively (see equations (TD] and [T2J), provided that the Mach number 
of the jet has values M w ~ 10 or larger. 



3. Simulations of the internal working surface of an HH jet 



Let us now consider the case of a jet model with a "top hat" density and velocity initial 
cross section, and an initial toroidal magnetic field cross section of the form 

r 



Br) 



Bn-. 



(15) 



This kind of magnetic field cross section has been used in many of the previous simulations 
of radiative MHD jets (see, e. g., iGardiner fc Frank! |2000| ). With this cross section for the 
jet beam, we have 

f£ i = -rj/r, (16) 



f 
■j i 



nr 
m 



-2r/rj 



and 



r-ni 
Jp 







rad 
V 



(17) 



In other words, the magnetic pressure force is directed towards the axis, and the gas pressure 
force (acting in the radial direction on the working surface jet material) is zero for the non- 
radiative case, and points outwards for the radiative case. 

We now compute models of a jet with this initial cross section, and an initial scale of the 
magnetic field Bq = (i. e., a purely hydrodynamic jet) and Bq = 5 /iG. The jet is injected 
with a constant density rij = 100 cm -3 , temperature Tj = 900 K and radius rj = 2 x 10 15 cm, 
and moves into a homogeneous, unmagnetized environment of density n env = 10 cm -3 and 
temperature T env = 9000 K. The injection velocity varies sinusoidally with time, with a 
period r = 20 yr, a half-amplitude of 150 km s _1 , and an average velocity of 300 km s _1 . 

For the two chosen values of the magnetic field (Bq = and 5 /iG, see above and 
Equation |T5"|) . we run both non-radiative simulations and simulations in which we include 
the coronal ionization equilibrium cooling function of iDalgarno fc McCrayl (119721 ). These 
simulations are run wit h the uniform grid, axisymmetric MHD code described in detail by 
De Colle fc Ragal (120061 ). The codes uses a second order up- wind scheme, which integrates 
the MHD equations using a Godunov method with a Riemann solver. The Riemann problem 
is solved using primitive v ariables and the magnetic field divergence is maintained close to 
zero using the CT method (lT6thll2000l ). The computational domain of (5, 1) x 10 16 cm (axial, 
radial) extent is resolved with 2000 x 400 grid points. A reflection condition is applied on 
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the jet axis and on the z = plane in the r > rj region. An outflow condition is applied in 
the remaining grid boundaries. 

The time-dependent ejection velocity of the jet leads to the formation of successive 
internal working surfaces that travel down the jet flow. It is possible to estimate the ratio 
(M/P) between the magnetic and pressure forces within the internal working surfaces by 
noting that the shock velocity (associated with the two working surface shocks) has a value 
v ~ 150 km s -1 . In other words, the value of the shock velocity is of th e order of the 



half-amplitude of the ejection velocity variability (see, e. g.. lRaga et al.lll990l ) 



With this value of v and the initial jet density and temperature, we can compute Ma = 
v/va ~ 31 (where va = 4.8 km s -1 for our B Q = 5 /xG value and our initial jet density), 
M w = 150 (for an assumed post-cooling sound speed of 10 km s" 1 ) and then we use equation 
( EED to obtain (M/P) nr ~ 0.02 and (M/P) rad ~ 10 3 . Therefore, the magnetic force should 
have little effect in the non-radiative simulations, and result in similar structures for the 
internal working surfaces in the cases of magnetized and non-magnetized jets. 

Figure 2 shows that our numerical simulations do show this effect. In this Figure, 
we show the density stratification obtained for non-radiative jets with Bq = (left) and 
B = 5 fiG (right) after a t = 90 yr time-integration. It is clear that though the details of 
the flow are affected by the presence of a toroidal magnetic field, the general features of the 
two working surfaces produced within the computational domain are quite similar in the 
magnetized and non-magnetized cases. 

The fact that (M/P) ra d ~ 700 (see above) implies that the magnetic pinch force should 
dominate the dynamics of the material within the internal working surfaces. Our two ra- 
diative numerical simulations (shown in Figure 3) do show this effect. In the magnetized 
simulation, the internal working surfaces become more strongly compressed towards the jet 
axis as they evolve (traveling away from the source), an effect not seen in the non-magnetized, 
radiative jet simulation (see Figure 3). 

Figure 4 shows zooms of the knot situated at z ~ 3.5 x 10 16 cm (the knot on the top 
half of the t = 90 yr time frames shown in Figures 1 and 2) for our four computed models. 
This Figure shows that in the non-radiative case, the two working surface shocks have a 
separation which is similar to the diameter of the jet, and that the density structures are 
very similar for the B = and 5 /iG models. 

As expected, much higher densities are obtained in the radiative jet simulations. In 
the radiative case, the working surface obtained from the Bq = 5 /iG model shows larger 
densities, a much higher concentration towards the jet axis and larger separations between 
the working surface shocks than the B = model. 
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More complex pro files for the magnetic fiel d and the pressure were explored in the past 



by several authors fe.g. iGardiner fc 



O 'Sullivan fc Ravi 12000 



Frank et al 



rank 



2000 



2000 



Gardiner et al. 2000; Stone fc Hardee 2000 



De Colle &: Ragal 120061 ). with results similar to 



the one obtained by the simple configurations of magnetic field and pressure presented here. 
As shown in Section 2.2, the expansion or contraction of the material in the working surface 
is nearly independent on the initial profile of the pressure. On the other side, the pre-shock 
magnetic field profile contributes to f r v ad (but not to and to the value of (M/P). For 

small values of the magnetic field (e.g. close to the jet axis) (M/P) < 1, and the gas pressure 
force dominates. 



Also, we note that IBegelmanl (119981 ) has studied the development of pinch instabilities 
in non radiative jets, as due to the presence of a toroidal magnetic field. He found that a 
condition necessary for the development of the pinch instability is 



dlnB <y8 

> 



dlnr 7/5 + 2' 



'11 



where (3 = MP/B 2 . In the case of a radiative working surface with a post-shock region with 
(3 <C 1 (corresponding to the condition (M/P) ra d 3> 1) this condition reduces to F m < 0. 



4. Conclusions 

It is a known result that internal working surfaces in radiative, MHD jets with a toroidal 
magnetic field configuration form dense, axial structures, which do not appear in unmagne- 
tized jets. We present a simple, analytic model with which we show that the strong jump 
conditions (applied to one of the working surface shocks) imply that the magnetic force 
dominates over the gas pressure force within a radiative working surface and that the gas 
pressure force is dominant for a non-radiative working surface (provided that one has a shock 
Mach number of at least M w ~ 10 and an Alvenic Mach number which does not exceed M 2 ). 

Interestingly, the radial dependence of the toroidal magnetic field within a radiative 
working surface depends only on the cross section of the pre-shock ram pressure p ra m( r ) = 
p(r)v 2 (r) impinging on the shocks. From equation dTUJ), we can see that if we have 

Pram (j ) 

that decreases towards the edge of the jet faster than 1/r 2 , the magnetic force within the 
working surface will be directed outwards, and will tend to increase the width of the working 
surface. 

We have run four simulations (with a top hat cross section for p ra m, that results in an 
axially directed magnetic pinch within the working surfaces), therefore in complete consis- 
tency with our analytic model. We find that in the non-radiative case the presence of a 
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toroidal magnetic field has very little effect on the structure of the internal working surfaces. 
We also find that for the radiative case, the presence of a toroidal magnetic field produces a 
strong axial compression of the material within the internal working surfaces (see Figure 4). 

The analytic model presented in this paper can then be used to decide what ram pressure 
and toroidal magnetic field cross section to use in a magnetized, radiative, variable jet 
simulation in order to produce internal working surfaces that show narrower or broader 
structures than what is obtained in non-magnetized jet simulations. This might be a valuable 
tool when trying to model the knots in specific HH jets, and might provide a possible method 
for constraining the strength and the configuration of magnetic fields within such objects. 
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Fig. I. — Schematic diagram showing an internal working surface produced by an ejection 
velocity variability (the source is to the left). In a frame of reference that moves with the 
velocity v ws of the working surface, material from the continuous jet segments enters the 
shocked layer from both the upstream and downstream directions. The B(r), P(r), v(r), 
p(r) radial cross section of the pre-shock region produces a B w (r), P w (r), p w (r) cross section 
within the shocked layer (in this shocked layer, the velocity along the jet axis is ~ in the 
reference frame moving with the working surface). The material in the working surface exits 
laterally, shocking against the jet cocoon. 




Fig. 2. — Density structures resulting from t = 90 yr time integrations of the non- magnetized 
(left) and magnetized (right) non-radiative jet models. The flow is injected from the bottom 
of the grid, and travels upwards, producing internal working surfaces (two of these are seen 
in the displayed time frames). The frames cover the full, (5,1) x 10 16 cm (axial, radial) 
computational domain. The densities are shown with the logarithmic gray scale given (in 
g cm~ 3 ) by the bar on the right. 
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Fig. 3. — Density structures resulting from t = 90 yr time integrations of the non- magnetized 
(left) and magnetized (right) radiative jet models. The frames cover the full, (5, 1) x 10 16 cm 
(axial, radial) computational domain. The densities are shown with the logarithmic gray 
scale given (in g cm -3 ) by the bar on the right. 
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Fig. 4. — Density stratifications of the knot seen in the upper half of the t = 90 yr time frames 
shown in Figures 2 and 3. The non-radiative models are shown on the top, and the radiative 
models on the bottom. The non- magnetized models are on the left, and the magnetized ones 
on the right. The displayed domain has an axial and radial size of 7.5 x 10 15 cm. The density 
of the non-radiative models is given (in g cm -3 ) by the bar on the top right, and the density 
of the radiative models by the bar on the bottom right. 



